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Equilibrium and absolute minimal states of Mumford-Shah 
functionals and brittle fracture propagation 

Marius Buliga* 
Abstract 



By a combination of geometrical and configurational analysis we study the 
properties of absolute minimal and equilibrium states of general Mumford-Shah 
functionals, with applications to models of quasistatic brittle fracture propagation. 
P I ■ The main results concern the mathematical relations between physical quantities 

as energy release rate and energy concentration for 3D cracks with complex shapes, 
seen as outer measures living on the crack edge. 

•*-> 
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1 Introduction 

> 

■ A new direction of research in brittle fracture mechanics begins with the article of 

Mumford & Shah [12] regarding the problem of image segmentation. This problem, 
which consists in finding the set of edges of a picture and constructing a smoothed 
version of that picture, it turns to be intimately related to the problem of brittle crack 
evolution. In the before mentioned article Mumford and Shah propose the following 
! variational approach to the problem of image segmentation: let g : Q C M 2 —> [0, 1] be 

the original picture, given as a distribution of grey levels (1 is white and is black), let 
u : Q — > R be the smoothed picture and K be the set of edges. K represents the set 
where u has jumps, i.e. u £ C 1 (Q\K,R). The pair formed by the smoothed picture u 
and the set of edges K minimizes then the functional: 

I(u,K) = [a \Vu\ 2 dx + [ (3 \ u-g \ 2 dx + . 
Jn Jn 

The parameter a controls the smoothness of the new picture u, (3 controls the 1? 
distance between the smoothed picture and the original one and 7 controls the total 
length of the edges given by this variational method. The authors remark that for /3 = 
the functional I might be useful for an energetic treatment of fracture mechanics. 

An energetic approach to fracture mechanics is naturally suited to explain brittle 
crack appearance under imposed boundary displacements. The idea is presented in the 
followings. 
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The state of a brittle body is described by a pair displacement-crack, (u, K) is such 
a pair if K is a crack — seen as a surface — which appears in the body and u is a 
displacement of the broken body under the imposed boundary displacement, i.e. u is 
continuous in the exterior of the surface K and u equals the imposed displacement uo 
on the exterior boundary of the body. 

Let us suppose that the total energy of the body is a Mumford-Shah functional of 
the form: 



The first term of the functional E represents the elastic energy of the body with the 
displacement u. The second term represents the energy consumed to produce the crack 
K in the body, with the boundary displacement uo as parameter. Then the crack that 
appears is supposed to be the second term of the pair (u, K) which minimizes the total 
energy E. 

After the rapid establishment of mathematical foundations, starting with De Giorgi, 
Ambrosio [8] , Ambrosio pQ , [2] , the development of such models continues with Franc- 



fort, Marigo [9], [TO], Mielke [II], Dal Maso, Francfort, Toader, [7], Buliga 0], [5], 



In this paper we introduce and study equilibrium and absolute minimal states of 
Mumford-Shah functionals, in relation with a general model of quasistatic brittle crack 
propagation. 

On the space of the states of a brittle body, which are admissible with respect to 
an imposed Dirichlet condition, we introduce a partial order relation. Namely the state 
(u, K) is "smaller than" (v, L) if L C K and E(u, K) < E(y, L). Equilibrium states for 
the Mumford-Shah energy E are then minimal elements of this partial order relation. 
Absolute minimal states are just minimizers of the energy E. 

Both equilibrium states and absolute minimal ones are good candidates for solutions 
of models for quasistatic brittle crack propagation. Usually such models, based on 
Mumford-Shah energies, take into consideration only absolute minimal states. However, 
it seems to me that equilibrium states are better, because it is physically sound to define 
a state of equilibrium (u, K) of a brittle body as one with the property that its total 
energy E(u, K) cannot be lowered by increasing the crack further. 

For this reason we study here properties of equilibrium and absolutely minimal 
states of general Mumford-Shah energies. This study culminates with an inequality 
between the energy release rate and elastic energy concentration, both defined as outer 
measures living on the edge of the crack. This result generalizes for tri-dimensional 
cracks with complex geometries what is known about brittle cracks with simple geom- 
etry in two dimensions. In the two dimensional case, for cracks with simple geometry, 
classical use of complex analysis lead us to an equality between the energy release rate 
and elastic energy concentration at the tip of the crack. We prove that for absolute 
minimal states (corresponding to cracks with complex geometry) such an equality still 
holds, but for general equilibrium states we only have an inequality. Roughly stated, 
such a difference in properties of equilibrium and absolute minimal states comes from 
the mathematical fact that the class of first variations around an equilibrium state is 
only a semigroups. 
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This research might be relevant for 3D brittle fracture criteria applied for cracks 
with complex geometries. Indeed, it is very difficult even to formulate 3D fracture 
criteria, because in three dimensions a crack of arbitrary shape does not have a finite 
number of "crack tips" (as in 2D classical theory), but an "edge" which is a collection 
of piecewise smooth curves in the 3D space. 

Aknowledgements. The author received partial support from the Romanian Min- 
istry of Education and Research, through the grant CEEX06-11-12/2006. 

2 Notations 

Partial derivatives of a function / with respect to coordinate Xj are denoted by fj . We 
use the convention of summation over the repeating indices. The open ball with center 
x G M n and radius r > is denoted by B(x, r). 

We assume that the body under study has an open, bounded, with locally Lipschitz 
boundary, reference configuration 0, C W 1 , with n = 1, 2 or 3. In the paper we shall use 
Hausdorff measures 7i k in W 1 . For example, if n = 3 then Ti n is the volume measure, 
H n is the area measure, Ti n ~ 2 is the length measure. If n = 2 then 7i n is the area 
measure, TL n ~ l is the length measure, TL n ~ 2 is the counting measure. 

Definition 2.1 A smooth diffeomorphism with compact support in £1 is a function 
(j) : — > £1 with the following properties: 

i) <j) is bijective; 

ii) (p and (f)~ l are C°° functions; 

Hi) (j) equals the identity map of Q near the boundary dfl: 

supp (idn — fa) CC H . 

The set of all diffeomorphisms with compact support in ft is denoted by D or T>{Vt). 

The set T){Q) it is obviously non void because it contains at least the identity map 
ida- Remark also that it is a group with respect to function composition. 

For any C°° vector field rj on £1 there is an unique associated one parameter flow, 
which is a function <p ■ I x & — > f2, where 7 C 1 is an open interval around G M, with 
the properties: 

fl) "it € I the function <p(t, •) = <f>t( m ) satisfies i) and ii) from definition 12.11 

f2) Vt, t' G J, if t — t' G / then we have (f>f o <f)^ 1 = cpt-t' > 

f3) V t G I we have 77 = fa o <f>^ 1 , where fa means the derivative oft>—>fa. 
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The vector field rj = generates the constant flow <f> t = idci- If r\ has compact 
support in then the associated flow 1 1— ► <pt is a curve in D. 

A crack set K is a piecewise Lipschitz surface with a boundary. This means that 
exists bi-Lipschitz functions (f a )aei-M, each of them defined over a relatively open 
subset D a of R+ _1 = {y G M"" 1 : y n _i > 0}, with ranges in M n , such that: 

K = U^ =1 f a (D a ) , 

if a + H then f a (D a \ OWf 1 ) n \ dRf 1 ) = . 

The edge of the crack K is defined by 

dK = U^/ a (D Q n dRIp 1 ) . 

We shall denote further by B r (dK) the tubular neighborhood of radius r of dfT, given 
by the formula: 

B r (dK) = U xedK B(x,r) . 

We denote by [/] = f + — f~ the jump of the function / over the surface K with respect 
to the field of normals n. 

3 Mumford-Shah type energies 

Definition 3.1 We describe the state of a brittle body by a pair (v,S*). The crack is 
seen as a piecewise Lipschitz surface S in the topological closure £1 of the reference 
configuration Q of the body and v represents the displacement of the body from the 
reference configuration. The displacement v has to be compatible with the crack , i.e. 
v has the regularity C l outside the surface S. 

The space of states of the brittle body with reference configuration Q is denoted by 
Stat(Q). 

The main hypothesis in models of brittle crack propagation based on Mumford-Shah 
type energies is the following. 

Brittle fracture hypothesis. The total energy of the body subject to the boundary 
displacement uo depends only on the state of the body (v, S) and it has the expression: 

E(v,S) = [ u)(Vv) dx + F(S;u ) . (3.0.1) 
Jn 

The first term of this functional is the elastic energy associated to the displacement v; 
the second term represents the energy needed to produce the crack S, with the boundary 
displacement uq as parameter. 

We suppose that the elastic energy potential w is a smooth, non negative function. 

The most simple form of the function F is the Griffith type energy: 

F(S; u ) = Const. ■ Area (S) , 

that is the energy consumed to create the crack S is proportional, through a material 
constant, to the area of S. 
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One may consider expressions of the surface energy F, different from (|3.0.ip . for 
example: 

F(v,5) = / 0(v+,v-,n) ds , 
Js 

where n is a field of normals over S , v + , v are the lateral limits of v on S with 
respect to directions n, respective — n and (j) has the property: 

<?Hv + ,v~,n) = 0(v~,v + ,-n) . 

The function eft, depending on the displacement of the "lips" of the crack, is a potential 
for surface forces acting on the crack. The expression (|3.0.ip does not lead to such 
forces. 

In general we shall suppose that the function F has the properties: 
hi) is sub-additive: for any two crack sets A, B we have 

F(AuB;u ) < F(A;u ) + F(B;u ) , 

h2) for any x € O and r > 0, let us denote by 5% the dilatation of center x and 
coefficient r: 

$r(y) =x + r(y-x) . 

Then, there is a constant C > 1 such that for any A C ft with F(A; Uo) < +oo 
we have: 

F(<£(A)nfi;uo) < Cr^FiA; u ) . 

The particular case F(A; uq) = GTC 71-1 (A) satisfies these two assumptions. In gen- 
eral these assumptions are satisfied for functions F(-;uq) which are measures absolutely 
continuous with respect to the area measure H n_1 . 

A weaker property than h2), is the property h3) below. We don't explain here why 
h3) is weaker than h2), but remark that h3) is satisfied by the same class of examples 
given for h2). 

For any icfi, let us denote by B(A,r) the tubular neighborhood of A: 

B(A,r) = U xeA B(x,r) . 

We shall suppose that F satisfies: 

h3) for any A C such that F(A; uq) < +oo, we have 

F(dB(A,r) n u ) 
limsup < +oo . 

r->o r 
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4 The space of admissible states of a brittle body 



Definition 4.1 The class of admissible states of a brittle body with respect to the crack 
F and with respect to the imposed displacement uq is defined as the collection of all 
states (v, S) such that 

(a) u = u on dQ \ S, 

(b) FCS U . 

This class of admissible states is denoted by Adm(F, Uq). 

An admissible displacement u is a function which has to be equal to the imposed 
displacement on the boundary of Q (condition (a)). Any such function u is reasonably 
smooth in the set f2 \ S u and the function u is allowed to have jumps along the set 
S. Physically the set represents the collection of all cracks in the body under the 
displacement u. The condition (b) tells us that the collection of all cracks associated 
to an admissible displacement u contains F, at least. 

For some states (u, S), the crack set S may have parts lying on the boundary of 0, 
that is S n <9f2 is a surface with positive area. In such cases we think about S D dQ as 
a region where the body has been detached from the machine which imposed upon the 
body the displacement uo. 

In a weak sense the whole space of states of a brittle body may be identified with 
the space of special functions with bounded deformation SBD(S7), see [3]. Indeed, to 
every displacement field u which is a special function with bounded deformation we 
associate the state of the brittle body described by (u, S u ), where generally for any set 
A we denote by A the topological closure of A. (Note that, technically, the crack set 
S u may not be a collection of surfaces with Lipschitz regularity.) 

On the space of states of a brittle body we introduce a partial order relation. The 
definition is connected to definition 14.11 and the brittle fracture hypothesis. 

Definition 4.2 Let (u, S),(v,L) £ Stat(£l) be two states of a brittle body with refer- 
ence configuration O. If 

(a) S CL, 

(b) u = v on dQ, \ L, 

(c) E(y,L)<E{vL,S), 

then we write (v, L) < (u, S). This is a partial order relation. 

There are many pairs (u, S), (v, L) G Stat(Q) such that (v, L) < (u, S) and (u, S) < 
(v,L), but u ^ v. Nevertheless such pairs have the same total energy E, the same 
crack set S = L, and u = v on d£l \ L. 

For a given boundary displacement uo and for given initial crack set K, on the set 
of admissible states Adm(uQ, K) we have the same partial order relation. 
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Definition 4.3 An element (u, S) G Adm(uQ, K) is minimal with respect to the partial 
order relation < if for any (v, L) G Adm(uQ, K) the relation (v, L) < (u, S) implies 
(eu,S)<(v,L). 

The set of equilibrium states with respect to given crack K and imposed boundary 
displacement vlq is denoted by Eq(uo,K) ant it consists of all minimal elements of 
Adm(uQ, K) with respect to the partial order relation <. 

An element (u, S) G Adm(uo, K) with the property that for any (v, L) G Adm(uo, K) 
we have E(u,S) < E(v,L), is called an absolute minimal state. The set of absolute 
minimal states is denoted by Absmin(uQ, K) . 

The physical interpretation of equilibrium states is the following. An equilibrium 
state (u, S) G Eq(iiQ,K) is one such that any other state (v, L) G Adm(uQ, K), which is 
comparable to (u, S) with respect to the relation <, has the property (u, S) < (v, L). In 
other words, equilibrium states are those with the property: the total energy E cannot 
be made smaller by prolongating the crack set S or by modifying the displacement u 
compatible with the crack set S and imposed boundary displacement uo. 

Absolute minimal states are just equilibrium states with minimal energy. 

Remark 4.4 There might exist several minimal elements of of Adm(uQ, K) , such that 
any two of them are not comparable with respect to the partial order relation <. 

For given expressions of the functions w and F, we formulate the following 

Equilibrium hypothesis (EH). For any piecewise C 1 imposed boundary displace- 
ment uo and any crack K the set of equilibrium states Eq(uQ,K) is not empty. 

Without supplementary hypothesis on the total energy E, the EH does not imply 
that the set of absolute minimal states Absmin(uQ, K) is non empty. Therefore the 
following hypothesis is stronger than EH. 

Strong equilibrium hypothesis (SEH). For any piecewise C 1 imposed boundary 
displacement uo and any crack K the set of equilibrium states Absmin(uQ, K) is not 
empty. 

5 Models of quasistatic evolution of brittle cracks 

We shall describe here two models of quasistatic brittle crack propagation, according 
to Francfort, Marigo [9], [TO], Mielke pi], section 7.6, or Buliga 0,0. At a first 
sight the models seem to be identical, but subtle differences exist. Further, instead 
of referring to a particular different model, we shall write about a general model of 
brittle crack propagation based on energy functionals, as if there is only one, general 
model, with different variants, according to the choice among axioms listed further. 
Whenever necessary, the exposition will contain variants of statements or assumptions 
which specializes the general model to one of the actual models in use. 

As an input of the model we have an initial crack set K C and a curve of imposed 
displacements t G [0, T] i— > Uo(i) on the boundary of $7, the initial configuration of the 
body. 
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We like to think about the configuration O as being an open, bounded subset of M n , 
n = 1,2,3, with sufficiently regular boundary (that is: piecewise Lipschitz boundary). 

The initial crack set K has the status of an initial condition. Thus, we suppose that 
d (R n \ $7) = d£l. For the same configuration Q we may consider any crack set K C Q 
as an initial crack. The crack set K may be empty. 

Remark 5.1 Models suitable for the evolution of brittle cracks under applied forces 
would be of great interest. Present formulations of the models of brittle crack propaga- 
tion allows only the introduction of conservative force fields, as it is done in ] 11^ or fWjj- 
The reason is that models based on energy minimization cannot deal with arbitrary force 
fields. In the case of a conservative force field it is enough to introduce the potential of 
the force field inside the expression of the total energy of the fractured body. Thus, in 
this particular case we do not have to change substantially the formulation of the model 
presented here, but only to slightly modify the expression of the energy functional. 

In order to simplify the model presented here, we suppose that no conservative force 
fields are imposed on f2 or parts of d£l. In the models described in [11] or [10J such 
forces may be imposed. 

Definition 5.2 A solution of the model is a curve of states of the brittle body t £ 
[0, T] i ^ (u(i)jjSt) such that: 

(Al) (initial condition) K C Sq, 

(A2) (boundary condition) for any t £ [0, T] we have u(t) = Uo(i) on d£l \ St, 

(A3) (quasistatic evolution) for any t £ [0, T] we have (u(t),St) 6 Eq(\io(t), St), 

(A4) (irreversible fracture process) for any t < t' we have St C St', 

(A5) (selection principle) for any t <t' and for any state (v, St) £ Adm(vLo(t') , St) we 
have E(y,S t ) > E(u(t'), S t >). 

From definition 14.31 we see that (A2) is just a part of (A3). The axiom (A2) is 
present in the previous definition only for expository reasons. 

The selection principle (A5) enforces the irreversible fracture process axiom (A4). 
Indeed, we may have severe non-uniqueness of solutions of the model. The axiom (A5) 
selects among all solutions satisfying (Al), (A4), the ones which are energetically 
economical. The crack set St does not grow too fast, according to (A5). For imposed 
displacement Uq(£ ), the body with crack set Sf is softer than the same body with the 
crack set St, for any t < t'. 

As presented in definition 15.21 the model has been proposed in Buliga [6]. 

In the models described in [TT], [9], [TO] we don't need the selection principle (A5) 
and the axiom (A3) takes the stronger form: 

(A3') (quasistatic evolution) for any t £ [0,T] we have (u(t),St) £ Absmin(uo(t) , St) ■ 



8 



6 The existence problem 



The existence of equilibrium, or absolutely minimal states clearly depends on the el- 
lipticity properties of the elastic energy potential w (as shown for example in [2], [3] 
or [U]). This is related to the existence of minimizers of the elastic energy functional, 
as shown by relation (jT.O.lj) further on. Some form of ellipticity of the function w is 
sufficient, but it is not clear if such conditions are also necessary. Much effort, especially 
of a mathematical nature, has been spent on this problem. 

In this paper we are not concerned with the existence problem, however. Our 
purpose is to find general properties of solutions of brittle fracture propagation models 
based on Mumford-Shah functionals. These properties do not depend on particular 
forms of the elastic energy potential w, but on the hypothesis made in the general model. 
As any other model, the one studied in this paper is better fitted to some physical 
situations than others. If some property of solutions of this model are incompatible 
with a particular physical case, then we must deduce that the model is not fitted for this 
particular case (meaning that at least one of the hypothesis of the model is not suitable 
to this physical case) . We are thus able to provide a complementary information to the 
one provided by the existence problem. See further the Conclusions section for more 
on the subject. 

7 Absolute minimal states versus equilibrium states 

The differences between the models come from the difference between equilibrium states 
and absolute minimal states. 

Absolute minimal states are equilibrium states, but not any equilibrium state is an 
absolute minimal state. 

Let us denote by (u, S) an equilibrium state of the body, with respect to the imposed 
displacement uo and initial crack set K. 

Consider first the class of all admissible pairs (v, S') such that S = S. We have, as 
an application of definition 14.31 then: 



Thus any equilibrium state minimizes the elastic energy functional (in the class of 
admissible pairs with the same associated crack set). A sufficient condition for the 
existence of such minimizers is the poly convexity of the elastic energy potential w. 

The elastic energy potential function w : M nxn (M) —* R associates to any strain 
F e M nxn (R) (here n = 2 or 3) the real value w(F) G R. If this function is smooth 
enough then we can define the (Cauchy) stress tensor as coming from the elastic energy 
potential: 



The variational inequality (|7.0.ip implies that in the sense of distributions we have: 





S C\n\K) . (7.0.1) 




div cr(u) = 
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and that on the crack set S we have 

<r(u) + n = <r(u)~n = , 

where the signs + and — denotes the lateral limits of cr(u) with respect to the field of 
normals n. 

7.1 Configurational relations for absolute minimal states 

We can also make smooth variations of the pair (u, S). Here appears the first difference 
between absolute minimal and equilibrium states. We suppose further that S \ K ^ 0, 
in fact we suppose that S \ K is a surface with positive area. 

If (v, L) E Adm(uQ, K) is an admissible state and cj) £ T> is a diffeomorphism of 
with compact support, such that K C (f>{K), then (v o cj) , 4>{S)) is admissible too. 

If (u, S) is an absolute minimal state then, as an application of definition 14.31 we 
have: 

E(u,S) < £(uof ^(S)) Vc/> G V , K C <f>{K) . (7.1.2) 

We may use (|7.1.2p in order to derive a first variation equality. 

We shall restrict further to the group T>(K) of diffeomorphisms (f> G T> such that 
supp {(p — id) nK = 0. Vector fields r\ which generate one-parameter flows in T)(K) are 
those with the property suppr] D K = 0. Further we shall work only with such vector 
fields. 

We shall admit further that for any smooth vector field rj there exist the derivatives 
at t = of the functions: 

t ^ I w^uo^ 1 )) dz , t ' ^ F(<k(K);ua) , 
Jn 

where 4>t is the one parameter flow generated by the vector field rj. The relation (I7.1.2P 
implies then: 

4, F(MS);uo) = - / W (V(uo^))d, . (7.1.3) 

at\t=o at\t=o Jq 

Let us compute the right hand side of (|7.1.3p . We have 



d_ 

dt\t. 



/ «)(V(uo cj) t 1 )) dx = / (-w(Vu) div r] + a(u) 1J (Vu) lfc (Vi|)fc j } dx . 
=o Jn Jo. 



For any vector field rj, let us define, for any x € S, X(x) = rj(x) ■ n(x), r] T (x) = 
rj(x) — A)(x)n(x), where n is a fixed field of normals over S. 

With these notations, and recalling that the divergence of the stress field equals 0, 
we have: 

-4- I wfVfuo^r 1 )) dx = [ HVujjAdH"- 1 + 
dt\t=o J n Yt Js L 

+ lim / {[w(Vu)]A - KuJ^Vu^} dH n ~ l . (7.1.4) 

r ^°JdB r (dS) 
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Definition 7.1 We introduce three kind of variations in terms of a vector field rj which 
generates an one parameter flow <f>t £ T){K): 

(a) (crack neutral variations) for r\ = on S; in this case we have 4>t(S) = S for 
any t, 

(b) (crack normal variations) for r) = An on S\K, with A:S-*Rd scalar, smooth 
function, such that \{x) = for any x £ K U dS, 

(c) (crack tangential variations) for r] ■ n = on S. 

For the case (a) of crack neutral variations the relation (|7.1.4|) gives no new information, 
when compared with (|7.0.ip . 

In the case (b) of crack normal variations, the relation (I7.1.4P implies 

| F(MK);u ) = I KVu^AdW"- 1 . 
dt \t=o Js 

In the particular case F(S;uq) = 'H n ~ 1 {S) we obtain: 

/ {KVu)] +H}\dH n - 1 = , 
Js 

where H = —div s n = — div n + rii jTiiiij is the mean curvature of the surface S. 
Therefore we have 

HVu)(z)] +H{x) =0 (7.1.5) 

for any x £ S \ K . 

In the case (c) of crack tangential variations, the relation (|7.1,4p implies 

| F(Msy,u ) = 

dt |t=o 

= Urn / {KVu)]A - ^uj^Vu)^]^^} dH 11 - 1 . (7.1.6) 

JdB r (dS) 

This last relation admits an well known interpretation, briefly explained in the next 
subsection. 

7.2 Absolute minimal states for n = 2 

Let us consider the case n = 2 and the function 

F(S;u ) = GH\S) , 

where TL l is the one-dimensional Hausdorff measure, i.e. the length measure. Let us 
suppose, for simplicity, that the initial crack set K is empty and the crack set S of 
the absolute minimal state (u, S) has only one edge, i.e. dS = {xq}. Let us choose 
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a vector field 77 with compact support in f2 such that r/ is tangent to S. The equality 
(|7.1.6p becomes then 

G r](xo) ■ t(x ) = lim f {[w(Vu)]r, • n - [a(u)„ (Vu)^]^-} dTT -1 , 

where t(x) is the unitary tangent m x £ K at K. If we suppose moreover that the 
crack S is straight near xq, and the material coordinates are chosen such that near xq 
we have r](x) = t(x) = (1,0), then the equality (|7.1.6p takes the form: 

G = lim / {HVu)]m - ^(u^Vu^K} dH^ 1 . (7.2.7) 

JdB r (x ) 

We recognize in the right term of (|7.2.7f) the integral J of Rice; therefore at the edge of 
the crack the integral J has to be equal to the constant G, interpreted as the constant 
of Griffith. 

The equality (|7.2.7p tells us that at the edge of a crack set belonging to an absolute 
minimal state the Griffith criterion is fulfilled with equality. 

7.3 Configurational inequalities 

For equilibrium states which are not absolute minimal states we obtain just an inequal- 
ity, instead of the equality from relation (|7.1.6p . Also, for such equilibrium states there 
is no relation like (|7.1.5|) between the mean curvature of the crack set and the jump of 
elastic energy potential. We explain this further. 

The reason lies in the fact that if (u, S) 6 Eq(uo,K) is an equilibrium state with 
S \ K having positive area, and (fi £ ~D(K) is a diffeomorphism preserving the initial 
crack set K, then we don't generally have the relation (I7.1.2p . 

Indeed, in order to be able to compare (u,5) with (u o </>(£)), we have to 
impose S C (fi(S). Only for these diffeomorphisms (fi £ *D{K) the relation (I7.1.2P is true. 
The class of these diffeomorphisms is not a group, like T>(K), but only a semigroup. 
Technically, this is the reason for having only an inequality replacing (I7.1.6p . and for 
the disappearance of relation (|7.1.5p . 

There is a necessary condition on the edge dS of the crack set S, in order to have a 
trivial vector field r/ which generates a one parameter flow (fit £ T^{K) with S C (fit(S) 
for any t £ [0, T] (with T > sufficiently small). This condition is dS \ K ^ 0. 

Thus, for (u, S) G Eq(iiQ, K) with S\K with positive area, and dS \ K / 0, we 
have 

E(u,S) < E(uo(fi~\(f) t (S)) VtG[0,T] , (7.3.8) 

for any one parameter flow (fit £ T^{K) with S C <fit{S) for any t G [0, T], 

In relation (j7.3.8H crack normal variations (case (b) of definition 17. ip are prohibited. 
But these type of variations led us to the relation (|7.1.5p . We deduce that for an 
equilibrium state (u, S) £ Eq(u.Q, K) , such that S\K has positive area, and dS\K ^ ill, 
the relation (|7.1,5p does not necessarily hold. 
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The crack tangential variations (case (c) of definition 17. f p are allowed in relation 
(|7.3.8p only for t > 0. That is why we get only a first variation inequality: 



- F(Msy,u ) > 

at \t=o 

> lim I {[w(Vu)]A - [a(u)y(Vu) ijfc ] %ni } dH"" 1 , (7.3.9) 

r ^°JdB r (dK) 

for any vector field r\ which generates one parameter flow 4>t £ F)(K) with 5 C <pt{S) 
for any t 6 [0, T]. 

The physical interpretation of relation (|7,3.9p is the following: the crack set S of 
an equilibrium state satisfies the Griffith criterion of fracture, but, in distinction with 
the case of an absolute minimal state, there is an inequality instead of the previous 
equality. We are aware of at least one example where this inequality is strict. This case 
concerns a crack set in 3D formed by a pair of intersecting, transversal planar cracks. 
Such a crack set has an edge (in form of a cross), but also a "tip" (at the intersection of 
the edges of the planar cracks. The physical implications of the inequality (|7.3.9p are 
that such a 3D crack may propagate in different ways, either along a crack tangential 
variation, or along a more topologically complex shape, by loosing its "tip". An article 
in preparation is dedicated to this subject. 

We may interpret the Griffith criterion of fracture, in the form given by relation 
(17.3.9p . as a first order stability condition for the crack S associated to the state of a 
brittle body. Surprisingly then, absolute minimal states are first order neutral (stable 
and unstable), even if globally stable (as global minima of the total energy). There 
might exist equilibrium states for which we have strict inequality in relation (|7.3.9I) . 
Such states are surely not absolute minimal, but they seem to be first order stable, if 
our interpretation of (17.3. 9p is physically sound. 



7.4 Concentration of energy from comparison with admissible states 

We can obtain energy concentration estimates from comparison of the energy of the 
equilibrium state (u, S) £ Eq(uQ,K) with other particular admissible pairs. 

Let xq £ Q be a fixed point and r > such that B(xo,r) C Q. We construct the 
following admissible pair (v r ,S r ): 



v r (x) 



u(x) if x G \ B{xq, r) 
if x £ Q n B(x ,r) 



S r = SUdB(x ,r) . 

We have then the inequality E(u, S) < E(v r , S r ), for any r > sufficiently small. We 
use the properties hi), h2) of F to deduce that for any xo G Cl and r > we have : 

f tu(Vu) dx < Cn n (x;u ) r n - 1 , (7.4.10) 

JB(x ,r) 

where £l n (xo; uq) is a number defined by 

f2 n (x ;u ) = F(dB(x , l);u ) . 
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In the case of Griffith type surface energy F(S;u ) = GH n 1 (S) we have 

Q n {xo; Uo) = GuJ n , 

with oj n the area of the boundary of the unit ball in n dimensions, that is oj\ = 2, 

U>2 = 27T, W3 = 47T 2 . 

This inequality lead us to the following energy concentration property for u: 

hmsup ; < Cn n (x Q ;u Q ) . 7.4.11 

The term from the left hand side of the relation (|7.4.1ip is the concentration factor 
of the elastic energy around the point xo. 

The relation (]7.4.1ip shows that the distribution of elastic energy of the body in 
the state (u, S) is what we expect it to be, from the physical viewpoint. Indeed, let 
us go back to the case n = 2. It is well known that in the case of linear elasticity in 
two dimensions, if (v, S) is a pair displacement-crack such that div <r(v) = outside 
S and <r(v) + n = <7(v)~n = on 5 then v behaves like ^fr near the edge of the 
crack, hence the elastic energy behaves like r _1 . We recover then the relation (|7.4.1ip 
for n = 2. 

The relation (|7.4.11|) does imply that elastic energy concentration has an upper 
bound, but it does not imply that the energy concentration is positive at the tip of 
the crack. In the case n = 2, for example, and for general form of the elastic energy 
density, the relation (|7.4.1ip tells us that if there is a concentration of energy (that 
is if the density of elastic energy goes to infinity around the point x in the reference 
configuration) then the elastic energy density behaves like r~ . But it might happen 
that the elastic energy density is nowhere infinite. In this case we simply have 

r Ibm w ( Vu ) dx n 

lim sup = = L) 

which is not in contradiction with (|7.4.11|) . 

From the hypothesis h3) upon the surface energy F we get a slightly different 
estimate. We need first a definition. 

Definition 7.2 For the equilibrium state (u, S) £ Eq(\iQ,K) and for any open set 
A C f2 we define: 



CE(u,S)(A) = lim sup 
r-»0 



J B ((dsnA,r)nn w ( Vu ) dx 



r.j^o ,. F(dB(dSnA,r);u ) 

CF(S;u )(A) = hmsup ■ 



^0 r 

The functions CE(u, S)(-), CF(S;uq)(-) are sub- additive functions which by well- 
known techniques induce outer measures over the a-algebra of borelian sets in O. 

The function CE(u, S)(-) is called the elastic energy concentration measure associ- 
ated to the equilibrium state (u,S). Likewise, the function CF(S;uq)(-) is called the 
surface energy concentration measure associated to (u, S) . 
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Theorem 7.3 Let (u, S) £ Eq(uQ,K) be an equilibrium state. Then for any open set 
A C we have 

CE(u,S)(A) <CF(S;u )(A) . 



Proof. We consider, for any closed subset A of the following admissible state 
(u r>A , S r ,A) given by: 



u r a{x) 



u(x) if x e fl\B(dSn A,r) 
if x £ fin B(dS H A, r) 



S rA = SUdB(dSDA,r) 



The state (u, S) is an equilibrium state and (u rjJ 4, S T; a) is a comparable state, therefore 
we obtain: 

f tu(Vu) dx < F(dB(dSn A,r);u ) . 

JB(dSnA,r)nn 

We get eventually: 

lB(dsnA,r)nn w ( Vu ) dx . ,. F{8B{dS n A r); u ) 
hm sup — < hm sup . U 

Theorem 17.31 shows that an equilibrium state satisfies a kind of Irwin type criterion. 
Indeed, Irwin criterion is formulated in terms of stress intensity factors. Closer inspec- 
tion reveals that really it is formulated in terms of elastic energy concentration factor, 
and that for special geometries of the crack set, and for linear elastic materials, we are 
able to compute the energy concentration factor as a combination of stress intensity 
factors. 



8 Energy release rate and energy concentration 

From relations (|7,1.3p , (|7.1.6p , we deduce that a good generalization of the J integral 
of Rice (which is classically a number) might a functional : 

d f 

rj , supp rj CC Q ^ / wiVfu.^V 1 )) dx , 

at \t=o J n 

where <pt is the flow generated by rj. 

Definition 8.1 For any equilibrium state (u, S) E -^(uq,^) and for any vector field 
r\ which generates a one parameter flow (fit £ T){K), such that (there is a T > with) 
S C <fit{S) for all t £ [0,T], we define the energy release rate along the vector field r] by: 

ER(u, S)(rj) = - ^ / »(V(uo^)) dx (8.0.1) 
at\t=0 Jq 
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Denote by V(K, S) the family of all vector fields rj generating a one parameter flow 
4>t £ V(K), such that there is a T > with S C 4> t {S) for all t £ [0,T]. Formally 
this set plays the role of the tangent space at the identity for the (infinite dimensional) 
semigroup of all <f> £ T>(K) such that S C (fr(S). 

Remark that ER(u, S)(r]) is a linear expression in the variable rj. Indeed, we have 

ER(u,S)(r]) = / {(T(Vu)ijUi t kr]k,j — tu(Vu) div 77} dx . 
Jn 

Nevertheless, the set V(K, S) is not a vector space (mainly because the class of all 
4> £ such that S C ^(S 1 ) is only a semigroup, and not a group). Therefore, the 

energy release rate is not a linear functional in a classical sense. 

Definition 8.2 With the notations from definition \8.1\ the total variation of the energy 
release rate in a open set D £\Q is defined by: 

I ER I (u,5)(D) = sup ER{u,S){ri) , (8.0.2) 

over all vector fields rj £ V(K,S), with support in D, supp n C D, such that for all 
x £ O we have \\rj(x)\\ < 1. 

The function \ ER \ (u, £)(•) is positive and sub-additive, therefore induces an outer 
measures over the a -algebra of borelian sets in 0. 

We call this function the energy release rate associated to (u, S) £ Eq(uo, K). 

The number | ER(u, S) \ (D) measures the maximal elastic energy release rate that 
can be obtained by propagating the crack set S inside the the set D, with sub- unitary 
speed, by preserving it's shape topologically. 

In the case n = 2, as explained in subsection 17.21 let x$ be the crack tip of the crack 
set S, and J the Rice integral. Then for an open set Dcfiwe have: 

- I ER(u, S) I (D) = J if the crack tip belongs to D, that is xq £ D, 

- I ER(u, S) I (D) = if the crack tip does not belong to D. 

For short, if we denote by 5xq the Dirac measure centered at the crack tip xq, we can 
write: 

I ER(u,S) |= J 5x . 

It is therefore the appropriate generalization of the Rice integral in three dimensions. 

Suppose that for any crack set L and boundary displacement uo the surface energy 
has the expression: 

F(S;u ) = GH n - l {S) . 

Then CF{S, mq){Q) is just G times the perimeter (length if n = 3) of the edge of the 
crack S which is not contained in K (technically, it is the Hausdorff measure 7i n ^ 2 of 
dS\K). 

There is a mathematical formula which expresses the perimeter of the edge of an 
arbitrary crack set L as an "area release rate". Indeed, it is well known that the 
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variation of the area of the crack set (f>t(L), along a one parameter flow generated by 
the vector field 77 £ V(K, L), has the expression: 

^ te- x {ih{S)) = I divtanv m n ~\x) , 

at |t=o J s 

where the operator divtan is the tangential divergence with respect to the surface S. If 
we denote by n the field of normals to the crack set S, then the expression of divtan 
operator is: 

divtanV = Vi,i ~ Vi,j n i n j ■ 

Further, the perimeter of dS \ K, the edge of the crack set S outside K, admits the 
following description, similar in principle to the expression of the elastic energy release 
rate given in definition 18.21 

H n - 2 {dS \K) = sup |y divtanV dH^ix) : 77 £ V{K, S), ViGl \\rj(x)\\ < 1 J . 

By putting together this expression of the perimeter, with relation (|7.1.6p . we obtain 
therefore the following proposition. 

Proposition 8.3 If for any crack set L we have F(L;u.q) = GTl n ~ 1 (L) then for any 
absolute minimal state (u, S) £ Absmin(\io, K) such that S\K 7^ we have 

I ER{u,S) I {fl) = CF{u,S)(n) . 



At this point let us remark that for a general equilibrium state in three dimensions 
(u,S) £ Eq(uQ,K) there is no obvious connection between the energy release rate 
I ER(u,S) I, as in definition 18.21 ant the elastic energy concentration CE(u,S), as in 
definition 17.21 

The following theorem gives a relation between these two quantities. 

Theorem 8.4 Let (u, S) £ Eq(uQ,K) be an equilibrium state of the brittle body with 
reference configuration 0,, and D <zVt an arbitrary open set. Then we have the following 
inequality: 

I ER(u, S) I (D) < CE(u,S)(D) . (8.0.3) 



Remark 8.5 For an arbitrary crack set L, we can't a priori deduce from the EH the 
existence of a displacement u' with (u ,L) £ Adm(uo, K) and such that for any other 
state (v, L) £ Adm(u.Q, K) we have 

/ tu(Vu') dx < kj(Vv) dx . 
Jn Jn 

From the mechanical point of view such an assumption is natural. There are mathemat- 
ical results which supports this hypothesis, but as far as I know, not with the regularity 
needed in this paper. Fortunately, we shall not need to make such an assumption in 
order to prove theorem \8.4\ 
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Proof. (First part) Let us consider an arbitrary vector field rj £ V(K, S), with com- 
pact support in D, such that for any x £ f2 we have ||?/(x)|| < 1. 
In order to prove the theorem it is enough to show that 

ER{u,S){ri) < CE(u,S){D) . (8.0.4) 

Indeed, suppose (|8.0.4p is true for any vector field ?? G V(K, S), with compact support 
in D, such that for any i£(!we have ||^(^)|| < 1. Then, by taking the supremum with 
respect to all such vector fields rj and using definition 18.21 we get the desired relation 
(l8TT3j) . 

The inequality f|8.0.4|) is a consequence of proposition 18. 61 and relation (|8.0.9[) . which 
are of independent interest. We shall resume the proof of theorem 18.41 by giving the 
proof of the inequality (|8.0.4|) . after we prove the before mentioned results. □ 

Let (ftt be the one parameter flow generated by the vector field rj. We can always 
find a curvilinear coordinate system (a±, a n -i, 7) in the open set D such that: 

- on the part of the edge dS n supp rj of the crack set 5 we have 7 = 0, 

- the surface 7 = t (constant) is the boundary of an open set Bt such that 

<f)t(S) \ S C B t C suppr] C D , 

- there exists T > such that for all t £ [0, T] we have 

B t C B(dSf]D,t)f]D , (8.0.5) 



where B(dS n D, t) is the tubular neighbourhood of dS n D, of radius t. 

Consider also the one parameter flow ipt, t > 0, which is equal to identity outside 
the open set D and, in curvilinear coordinates just introduced, it has the expression 

V>t(x(a;,7)) = x(ai,t + 7) . 

Notice that ipt(ty = ^ \ -Bt- We shall use these notations for proving that the elastic 
energy concentration is a kind of energy release rate, after the following result. 

Proposition 8.6 With the notations made before, we have: 

lim-f/ w(Vu)dx-[ w(V{u o tp' 1 )) dx ) =0 . (8.0.6) 
t \Jn\Bt Jn\B t J 
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Proof. Recalling that ipt(£l) = f2 \ B t , we use the change of variables x = tpt{y) to 
prove that (|8.U.6|) is equivalent with 

lirn^ U (^(Vu(y)(V^) _1 (y)) - w((Vu)(My))) detVMv) dy\ =0 . 

The previous relation is just 

^ / ( W (Vu( 2 /)(V^)- 1 (2/))-^((Vu)(Vt(y)))detVV t (y) dy =0 . (8.0.7) 
at |t=o Jq 

We shall prove this from (u, S) G £ , g(uo,-ftr) and from an approximation argument. 
Notations from subsection 17.11 will be in use. 

Denote by u the vector field which generates the one parameter flow tpf Let us 
compute, using integration by parts: 

^ / (w(Vu(y)(Vtl>tr\y))-w((Vu)(i> t (y))))detV^ t (y) dy = 
at \t=o J n 

= / {vijU-ijkUk + °ijUi,kUk,j) dy . (8.0.8) 
Jn 

For any 7 > 0, sufficiently small, choose a smooth scalar function / 7 : f2 — > [0, 1], 
such that: 

(a) f{x) = for all x G 5 7 , / 7 (x) = 1 for all a? G O \ B 2l , 

(b) as 7 goes to we have: 

lim / p (<7ijUjjfcu; fe + (TijUi^Ukj) dy = (crijUi t j k u k + ^ij^k^kj) dy , 

lim / f]aiju ik uj k dy = . 
For all sufficiently small 7 > it is true that: 

J (vijUijkuJ + aiju i>k u2^ dy = 

Thus, from (a), (b) above we get the equality: 

lim / (aijUukuZ + aijUi kU?-) dy = ((TijUu k u k + cr^U; Mk,j) dy • 
~f->0Jn \ ,J ' Jn 

Recall that (u, S) is an equilibrium state, therefore the stress field a = cr(Vu) has 
divergence equal to 0. Integration by parts shows that for any sufficiently small 7 > 
we have: 

J (°ij u i,jkul + (TijU i>k ul^ dy = J -a ijd (u ijfc w 7 ) dy = . 
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We obtained therefore the relation: 

/ ( a ij u ijfc^fe + &ijUi,kVk,j) dy =0 . 
Jn 

This is equivalent to relation (18.0. 7ft . by computation ()8.Q.8[) . □ 

A straightforward consequence of (|8,0.6p is that the elastic energy concentration is 
related to a kind of configurational energy release rate. Namely, we see that 

limsup- / w(Vu) dx = 

= limsup- ( / w(Vu)dx- [ u;(V(u o ^T 1 )) dx ) . (8.0.9) 
t yj n Jn\B t J 

We turn back to the proof of theorem 18.0.31 Recall that what it is left to prove is 
relation ([QTIjl . 

Proof of (I8.0.4|) . By construction, for all sufficiently small t > we have: 

If If 

- / w(Vu) dx > - Hf(Vu) dx . 

t JB(dS,t)nD t J Bt 



because B t C B(dS, t) D D. We write the right hand side member of this inequality as 
a sum of three terms: 

i r 

dx 

'B t 



7 / «>(Vu) 

* JB t 

j(J w(Vu)dx-J w(V(uo (f)- 1 )) dx^j + 



+ -( / w(V(uo^)) dx - [ lufV^o^ 1 ))^ + 
* \Jn Jn\B t J 

+ -([ w(V(uoV t -1 )) dx- I io(V(u)) dx) . 
1 \Jn\Bt Jn\B t ) 



in\B t Jn\B t 

As t goes to 0, the first term converges to EC(u, S)(rj) and the third term converges 
to by proposition 18.61 We want to show that 

lim- | / w(V(uo^)) dx- / w(V(u o t/^ 1 )) dx ) =0 . (8.0.10) 
t \Jn Jn\B t J 

The proof of this limit is identical with the proof of proposition 18.61 Indeed, in that 
proof we worked with the one parameter flow tpt generated by the vector field u. This 
one parameter flow is a semigroup (with respect to composition of functions), but after 
inspection of the proof it can be seen that we only used the following: for any x £ £1\S 

lim^t(x) = x and — ipt(x) = u(x) 

t >0 Clt |i^0 
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Therefore we can modify the proof of proposition 18.61 by considering, instead of rfe , the 
diffeomorphisms \t defined by: 

At = ipt ° (Pt 1 ■ 

The rest of the proof goes exactly as before, thus leading us to relation (|8.0. lQj) . 
Eventually, we have: 



CE(u,S){D) = limsup- [ mj(Vu) dx > 

t->o t J B(dS,t)r\D 

> limsup- / io(Vu) dx = ES(u, S)(rj) + 

t^O t JB(dS,t)C\D 
+ lim-( I w(V(u o dx - I w^uo^ 1 )) di + 

+ lim-( / wj(V(uo Vr 1 )) dx ~ I u>(V(u))drj = ES(u,S)(rj) 
t \Jn\Bt Jn\B t J 

and (|8.0.4p is therefore proven. □ 



9 A constraint on some minimal solutions 

Let us consider now a solution of the model of brittle crack propagation described in 
section [5j More precisely, for given boundary conditions uo(t) and initial crack set K, 
we shall call a solution (u(t),St) G Eq(\io(t), St) of the model described by axioms 
(Al),..., (A5), by the name "equilibrium solution". Likewise, a solution (u(t), <St) G 
Absmin(vLo(t) , St) of the model described by axioms (A1),(A2),(A3'),(A4), will be 
called a " minimal solution" . 

We shall deal with a minimal solution (u(t),St) € Absmin(uo(t) , St) for which the 
crack set St propagates smoothly, without topological changes. Namely we shall suppose 
that there exists a vector field rj with compact support in $7, such that for all t £ [0, T] 
we have St = <f>t(K), where <pt is the one parameter flow generated by r/. 

Because the problem is quasistatic, time enters only as a parameter, therefore we 
may suppose moreover that for all x £ we have rj(x) < 1. 

At each moment t G [0, T] we shall have r] o (j) t £ V(K, St)- 

Theorem 9.1 Suppose that for any crack set L and boundary displacement uq the 
surface energy has the expression: 

F(S;u ) = GH n - l {S) . 

Let (u(t), St) £ Absmin(uo(t), St) be a minimal solution, with So = K, such that exists 
a vector field n with \\r](x)\\ < 1 for all x £ Q and for all t £ [0, T] we have St = <pt{K), 
where 4>t is the one parameter flow generated by rj. 
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Then for any t £ [0, T] and any open set D C we have the equalities: 



ER(u(t),MS)) I (D) 



EC(u(t),MS))(D) 



CF(MS);u (t))(D) 



GH n - 2 {dS\K) 



(9.0.1) 



Proof. Theorems l8.4l and l7.3l tell us that for any open set D C f2, and for any t £ [0, T] 
we have 



We deduce that for any open set D C O, and for any t € [0, T] the string of equalities 
(19^011) is true. □ 

This result is natural in two dimensional linear elasticity. Nevertheless, in the 
case of three dimensional elasticity, the constraints on the elastic energy concentration 
provided by theorem 19.11 might be too hard to satisfy. 

Indeed, from (|9,0.ip we deduce that in particular the elastic energy concentration 
has to be absolutely continuous with respect to the perimeter measure of the edge of 
the crack. 

10 Conclusions 

We have proposed a general model of brittle crack propagation based on Mumford-Shah 
functionals. We have defined equilibrium and absolute minimal solutions of the model. 

By a combination of analytical and configurational analysis, we defined measures 
of energy release rate and energy concentrations for equilibrium and absolute minimal 
solutions and we have shown that there is a difference between such solutions, as shown 
mainly by theorems 17.31 18.41 and 19.11 
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